Electronic zero modes of vortices in Hall states of gapped graphene 

Gordon W. Scmcnoff 

Department of Physics and Astronomy, University of British Columbia 

6224 Agricultural Road, Vancouver, British Columbia V6T IZl 

Recent observation of a metal-insulator phase transition in the u — Hall state of graphene has 
inspired the idea that charge carriers in the metallic state could be fractionally charged vortices. 
We examine the question of whether vortices in particular gapped states of graphene and subject 
to external magnetic and pseudo-magnetic fields could have the mid-gap zero mode electron states 
which would allow them to be charged. 
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1. INTRODUCTION 

It is now established that, at sufficiently low temper- 
ature and strong magnetic field, the four-fold degener- 
ate Landau level which resides at the apex of the Dirac 
cone in graphene is split into four sublevels [i|-0. At 
least some of this splitting is attributed to spontaneous 
breaking of the sublatticc symmetry and generation of a 
mass gap in the electron spectrum. The further recent 
observation of Kosterlitz-Thouless scaling of the resistiv- 
ity near the metal-insulator transition which occurs as 
the strength of the field is varied [a]-[3 has inspired a 
number of possible explanations, perhaps the most inter- 
esting of which has fractionally charged vortices as the 
carriers of electric charge in the metallic phasej9||10|[ll|. 
This scenario would give another concrete physical ex- 
ample beyond polyacetylene [I4] of topological solitons 
with fractional quantum numbers. 

The vortices arc topological defects in the condensate 
which creates the mass gap. In a system with particle- 
hole symmetry, electron energy states above and below 
the middle of the mass gap are paired and fractional 
charge can only arise from unpaired mid-gap states [ij]. 
The possibility that the interaction of vortices with rel- 
ativistic fermions can result in mid-gap states has been 
known for a long time [ij]-[l3]- The novel feature of 
the Hall system is the presence of a macroscopic mag- 
netic field where the electron spectrum would normally 
be concentrated in Landau levels. In the absence of the 
mass condensate, or other symmetry breaking, a four-fold 
degenerate (two spins and two v alley s) Landau level re- 
sides precisely at the Dirac point [1^ . With a parity and 
time-reversal invariant mass gap 2m, this Landau level is 
split into two levels which are displaced to the positive 
and negative mass thresholds, E = ±m [l3| and, with- 
out defects, the mass gap contains no electronic states 
at all. It is then reasonable to ask whether a vortex de- 
fect in the condensate which creates the mass gap can 
bind a mid-gap electronic state. This question has been 
studied recently by Hou ct.al. [l^l and Herbut [11| who 
found, surprisingly, that the answer is yes. In fact, the 
existence of mid-gap states attached to vortices is rela- 
tively insensitive to the presence of a magnetic field. In 
the following, we shall examine this phenomenon in more 
detail. We will show that the fermion-vortex system in 
simultaneous magnetic and pseudo-magnetic fields gen- 



erally has \n\ zero modes, where n is the vorticity, even 
when the magnetic field is a constant so that its flux di- 
verges with the volume of the system. The exception is 
where the pseudo-magnetic flux inside a disc of large ra- 
dius r — >■ 00 grows faster than r and, in addition, when it 
is larger than the magnetic flux at large r. In that case, 
the number of zero modes no longer depends on the vor- 
ticity, instead there are two infinitely degenerate Landau 
levels at zero energy. 

We shall also show that, when the pseudo-magnetic 
field falls off at least as fast at l/r^ at r — !• 00, an exter- 
nal magnetic field has no effect on the 77-invariant of the 
graphene Dirac operator perturbed by a constant mass 
term which upsets the particle-hole symmetry. 77 is a 
topological invariant and it is of physical interest because 
it is proportional to the electric charge of a system of 
electrons which is governed by the single-particle Hamil- 
tonian in question and in a state where all of the negative 
energy levels are filled and all of the positive energy lev- 
els are empty. This result indicates that the fractional 
charge of a vortex is indeed a topological invariant and 
that it is independent of the external magnetic field. 

Finally, we shall find interesting correlations between 
the sublatticc components of the electron mid-gap state 
wavefunctions and its internal quantum numbers, such 
as spin or valley assignment. These correlations are 
analogous to the Dirac point Landau level in ungapped 
graphene, where, for one valley the Landau level zero 
mode wavefunctions all live on one of the two sublattices 
and the other valley zero mode wavefunctions live on the 
other sublatticc, which sublatticc depends on the sign of 
the external magnetic field. Zero modes in a vortex back- 
ground also follow a pattern that is similar to this. We 
will elaborate in later Sections. 

A number of scenarios for mass gap generation have 
been proposed [13|-[33l. To support topological vortices, 
the mass condensate should have an S0(2) symmetry. 
Two proposals have appeared in recent literature. One 
is antiferromagnetic order with an easy plane [3J]. For- 
mation of antiferromagnetic order would be a result of 
the short-ranged on-site Coulomb repulsion of electrons 
whose role is enhanced in a strong magnetic field [231 • The 
Zeeman coupling of electron spins to the magnetic field 
would result in an "easy plane" where the antiferromag- 
netic order parameter lies in the plane that is perpen- 
dicular to the applied field. This is easy to demonstrate 



using second order degenerate perturbation theory with 
the Zeeman coupling as a perturbation and the possi- 
ble orientations of the antiferromagnetic condensate as 
the set of degenerate ground states. The other scenario 
with SO (2) symmetry is a bond density wave called the 
Kekule distortion which is a result of the electron-phonon 
coupling |26l| . jOJ j34l| - |36l| . It was shown to have vortex de- 
fects which, when coupled to the graphene electron, bind 
an unpaired mid-gap state. 



2. GRAPHENE HAMILTONIAN WITH A MASS 

GAP 

We shall begin by considering graphene at energies 
near the Dirac points, and with a mass generated by anti- 
ferromagnetic order. The bipartite honeycomb lattice of 
graphene is a superposition of two triangular sublattices. 
Antiferromagnetic order gives each of the two spin states 
of the electron a charge density wave - spin up electrons 
have higher density on one of the two sublattices and 
spin down electrons on the other sublattice. This breaks 
the symmetry under interchanging the two sublattices 
and the charge density wave generates a parity-invariant 
mass term in graphene. [l8| . The electrons are described 
by the Hamiltonian 
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thvpff ■ D -i 
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(1) 
This Hamiltonian is an 8 x 8 matrix. The 2x2 structure 

which is displayed explicitly refers to the two graphene 
valleys. The four additional dimensions which are sup- 
pressed are a direct product of two dimensional matrices: 
one on which is described by the pseudo-spin Pauli ma- 
trices ct" act and one by the spin Pauli matrices r". An 
alternative way of presenting this Hamiltonian is by writ- 
ing the 8x8 matrices as ordered direct products of the 
three 2x2 matrices with unit and Pauli matrices for 
pseudo-spin (I, ct"), spin (X, r'') and valleys (1,7]'^), in 
that order. In terms of these matrices, ([1]) has the form 
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Here, the covariant derivative is D 



the electromagnetic vector potential and vp w -^ 



i-A with A 

c 

is the 
graphene Fermi velocity. The magnetic field is Ba = 

V X i. 

It is useful to remind the reader that, the upper com- 
ponent of the spinors on which the pseudo-spin matrices 
(X, (t") act has support on one of the graphene sublat- 
tices, the lower component on the other. 

In ([1]) and , the masses of each spin polarization are 
±|77i|, the eigenvalues of TOaT°. If the antiferromagnetism 
has an "easy plane" , only two of the the three spin ma- 
trices can be used to describe possible orientations of the 
condensate, for example, rh ■ t = rriiT^ + TO2T^. We will 



be interested in the circumstance where nia can depend 
on spatial position and can have a vortex profile, with 
asymptotic limit 

lim [mi (r, 6) + im2 (r, 6)]^ me*"^ + 01-) (3) 

r— >-oo \ ^ / 

with (r, 6) the polar coordinates of the plane and m a 
constant. 

If the mass term in ([1]) were absent, the existence of two 
valleys and two spin states give the Hamiltonian an SU(4) 
symmetry. The mass term reduces this to an SU(2) sym- 
metry which mixes the valleys. [37| It is generated by 

I (8) X ® 77^ , cr^ ® r^ (g) 77^ cr^ «) T^ (8) 77^ (4) 

Before we proceed, we observe that, we can redefine 
the matrix components of the Hamiltonian ^ by conju- 
gating it with the matrix 



I(E}I®-{I + r]^)+a^(l)T^(E} -{T-if) 



then interchanging the spin and valley labels, r" o 77" 
and then conjugating with 

the resulting Hamiltonian is identical to the one with a 
mass term arising from the Kekule distortion that was 
discussed in |2a and 3J]. This replacement amounts to 



an SU(4) rotation and other orientations of the SU(4) 
breaking terms in the Hamiltonian are possible as well. 
We also recall that in Ref . [3J] , when the mass is oriented 
as coming from a Kekule distortion, they add a pseudo- 
magnetic gauge field which couples to the Kekule conden- 
sate. In the present case where the condensate is antifer- 
romagnetic order, the analogous gauge field would modify 
the covariant derivative as V — i-A — ;> V — i-A — iVr^ 

c c 

and it would add a term 



' ®rf Vi+a'^ 



'Xy, 



in ([2]). In the following, we shall include this gauge field 
in our analysis. In the case where the condensate is anti- 
ferromagnetic, it amounts to gauging the symmetry that 
rotates the antiferromagnetic order parameter, as might 
be done to describe a spin liquid. 



3. ZERO MODES 

Since the valley degrees of freedom in ^ are de- 
coupled, we can perform our analysis for each valley 
separately. Consider the Hamiltonian describing elec- 
trons in one of the valleys (for simplicity, in units where 

h^ia' (g) XVj + cr* ® lAi + (j"- (g) T^Vi + (T^ (E) t" rua (5) 



As we have suppressed the valley degree of freedom (77" 

are missing), this is a 4 x 4 matrix differential operator. 

We will search for solutions of hipQ — 0. We note that 

the Hamiltonian ([S]) anticommutes with the matrix a^ Cg) 







(6) 



Zero modes of h can also be chosen to be eigenvec- 
tors of (j^ (g) T^ and we shall show below that all of 
the zero modes are indeed eigenvectors with eigenvalues 
cr'^ (g) r"^ = — sign(77i). This implies that the sublattice 
position and the spin orientation of the zero mode are 
correlated. One sublattice has a^ = 1 and on that sub- 
lattice the spin is polarized as t'^ = sign(n). The other 
sublattice has cr'^ — — 1 and the opposite spin polariza- 
tion T^ = — sign(7i). Furthermore, we shall find that, 
close to the vortex core, the zero modes have support on 
both sublattices. On the other hand, far from the vor- 
tex core, the zero modes have support on only one of the 
sublattices and therefore have one particular orientation 
of the spin. Which of the two possible sublattices is cho- 
sen depends on the orientation of the external magnetic 
field. Also, we note that the spin polarizations that we 
are discussing here are always orthogonal to the spin po- 
larization of the condensate which is in the r^-r^-plane. 

Also note that, if we include the second valley, its 
Hamiltonian differs from ([S]) by conjugation with the ma- 
trix I^a'^(E)T^ which interchanges the sublattices. Thus, 
for the other valley, there are also \n\ zero modes. They 
also have a peculiar spin polarization to sublattice rela- 
tionship which is precisely the opposite of the one in the 
first valley. 

We will begin with the special case of a rotation- 
ally covariant vortex with vorticity n: (7711,7712) = 
m{r){cos n9, sin n9) where m{r) ^ 7?7 at r — >■ 00 as 
in ^ and 777,(0) = 0. We shall also assume that the 
magnetic fields Ba and By = V x F depend only on 
r. We will use the gauge Ar — and Vr = where 
Ae{r) = J^r'dr'BA{r') and Ve{r) = J^ r'dr'Byir'). For 
example, if the magnetic fields are constants, the vector 
potentials would be Ae{r) = ^r^ and Ve{r) = ^r'^. 
When written in polar coordinates, the Dirac equation 
for a zero mode is 



m[r)e 



T u + le [Or oe — 



Ag + T^Vg 



V = 0(7) 



le'" I dr H — d( 
r 



u-m{r)e'^^ \'v = 0{%) 



The spinor components u and v are eigenstates of pseudo- 
spin a^ ® I with eigenvalues -1-1 and —1, respectively. 
Formally, they can be defined as 



u=\{X + a') 



)I t/iq 



1 



)I Ipo 



These project two two-component spinors, u and v, from 
the four-component spinor t/jq. u and v each have two 



spin components on which the Pauli matrices t" act. We 
make the Ansatze 



u{r, 6) = e*^« r^ u{r) 
v{r,e) = e-'''%{r) 

Plugging this into ^ and dH) yields 



(9) 
(10) 



r(r)e'(^+k+nr^+^Wu + ifo^ 



k + Ae+T^Ve 



w = 



,e^ii+k+nr-^+i}e f dr - L^LtlEll] u - m{r)i = 



Now, we shall have to choose the spinors u and v to 
be eigenvectors of t^, both with the same eigenvalue 
-sign(n), 

T u = — sigma(r7)M , t v ~ — sign(n)'D 

(The reason for this choice of sign will become clear when 
we find the possible values of k and £ which both will have 
to be positive and have to obey (jlip below.) u and v now 
ech have only one non-zero component. 
Then, the 6'-dependent phases cancel if 



£ + k^ \n\ -1 



fll) 



which we shall assume holds. Consistency of (fTTj) is a 
result of the correct choice of sign of the eigenvalues of 

What remains is 

f \- , -ffi k + Ag - sign{m)Ve\ . . 

m[r)u + I [ Or \ V = (12) 



,[dr-^—^l^^^^^^^^!^]u-m{r)i = (13) 



Assuming that Ag, Vg and 777(r) go to zero at small 
r, by examining (fT2|) and (fT3ll at r ~ 0, we see that the 
small r behavior of the wave-functions must be 



u{r) 
v(r) 






(14) 
(15) 



We shall need both of these solutions in order to match 
the solution at r — > 00 where there will generally only be 
one solution. Both solutions are normalizable at r = if 
i? = 0, 1, 2, ... and if A; = 0, 1, 2, .... However, remembering 
that k and i are related by (|11|) . we see that the allowed 
solutions are 

(£,fc) = (0,|77|-l),(l,|r7|-2),...,(|77|-l,0) , |r7|>l 

(16) 
We conclude that there are exactly |n| possible solutions 
where both behaviors are normalizable at the origin. 

The spinor components u and 7; are then both eigen- 
vectors of spinZcgJT^ with the same eigenvalue: I®t^u = 
— sign(ri)7i and I ® t^v = — sign(7i)'D. Recalling that u 
and u differ by a factor of Z ® r^ which flips the sign of 



the eigenvalue oi T ^ t^ , we conclude that u and v are 
eigenvectors of spin with 



T®t'^u = +sign(ri)w 
I®T^v — — sign(n)i' 



(17) 
(18) 



Also, we recall that u and v are eigenvectors of pseudo- 
spin a^ ® X: 



(7^ (g}Iv 



u 

V 



(19) 
(20) 



Combining the eigenvalues tells us that u and v are eigen- 
vectors of a^ (E) T^ with the same eigenvalue, sign(n): 



T u — sign(7n)u 
) TV = sign(?7i)w 



or, more succinctly 



r^Vo 



sign(TO)V'o 



(21) 
(22) 



(23) 



Remember that a^ (E) t^ anticommutes with the Hamil- 
tonian, so we expect that the zero mode wave-functions 
ipo are are also eigenvectors of this matrix. Here, we see 
explicitly that is is the case and the eigenvalue is deter- 
mined by the sign of the vorticity. This establishes the 
basis for the discussion after equation ^ above. 

Now we shall consider the large r limit of ([T2|) and (fT2|l . 
If Ag and Vg grow slower than r at r -> oo (meaning that 
Ba and By fall off faster than 1/r), the solution at large 
r is identical to the one found by Jackiw and Rossi [Ij] 
and their conclusion that there are exactly \n\ zero modes 
applies here. To generalize this, we will assume that at 
least one of Ag and Vg grow faster than r as r — >■ cx). 
(This means that BA{r) or Bv{r) fall off slower than 
1/r and one of both of the fluxes J (PxBa or J d^xBv 
diverge.) Then (fT^ and (fT^ are solved at large r by 



u{r) 
v{r) 



J; '^iAeir')-signi,n}Veir'}) ~ 



Wo 



Jo ^(-Ae(r')-sign(Tn)Ve(r')) ~ 



VO 



(24) 
(25) 



In addition, in the large r regime, the solutions re- 
side on either one or the other sublattice (either u oi v 
is non-zero), depending on the sign of Ag ( and there- 
fore the magnetic flux (pA = ^ J (PxBa- Since the spin 
and the pseudo-spin are correlated, all zero modes in the 
asymptotic region have one particular spin polarization - 
the eigenvalue of r^ is — sign((/)yi)sign(n) - the zero mode 
in the asymptotic region is polarized in one particular di- 
rection, orthogonal to the easy plane of the condensate. 
The direction depends on the sign of the total magnetic 
flux and the sign of the vorticity. 

Now, consider the case where \Ag\ < \Vg\. To under- 
stand this case more clearly, we observe that, if V{x) is 
in the Coulomb gauge, V • F = 0, we can rewrite h as 



5t' XV {x) 



{^ 



)ZV, 



+a^ 



r'^m^e^' 



E)IAi 

r^Xvix) 



(26) 



where Xv{x) = ^ J cPx' In \x — x'\By{x'). The zero 
mode problem for h can be solved by flnding potential 
zero modes when V ~ and testing their normalizability 
at r = oo when they are multiplied by the additional fac- 
tor e^'^ ^"^ xvix)^ jjj ^jjg gg^gg lY^g^i -^ve are interested in, 
when e -*^^ dominates the asymptotics, there are always 
two normalizable solutions with correlated pseudo-spin 
and spin so that a^ (8) t^xv{oo) > 0. When we begin 
at small r with one of the spinors with behavior in (J14p 



or (115 



u(r) 











and extrapolate 



to large r, they should become linear combinations of 
the two normalizable solutions there. In this case, there 
are an infinite number of zero modes for £ = 0, 1,2, ... 
or k = 0,1,2,... which comprise a two- fold degenerate 
Landau level at zero energy. We conclude that, when 
\Ae\ < \Ve\, there is a Landau level at the Dirac point. 
Here, by "Landau level" we mean an infinite set of de- 
generate states and they exist whenever Bv{r) falls off 
slower than 1/r at r — > oo. 



4. INDEX AND ??-INVARIANT 



We must now choose mq and vq so that the wavefunctions 
are normalizable. 

First consider the case where \Ag\ > \Vg\. Then, if 
Ag > ai large r, we set uq ^ and (P5|) is normalizable 
at r — ?► oo. On the other hand, if Ag < 0, we set vq = 
and (p4|) is normalizable. Thus we see that, in the case 
where \Ag\ > \Vg\, there is always one normalizable solu- 
tion of the Dirac equation at large r. Once it is given the 
appropriate angle dependence, e^^^T^u{r) and e~''^^z;(r) 
with one of the \n\ allowed values of k and £, and con- 
tinued to small r, it will become a linear combination of 
the two normalizable solutions that we found at r ^ 0. 
We conclude that, when \Ag\ > \Vg\, there are exactly 
\n\ zero modes, independent of the specific asymptotic 
behavior of Ag. This is in line with the conclusions in 
Refs. [13 and [n I . 



The index of h is the number of zero modes which 
have eigenvalue of cr^ (g) t'^ = +1 minus the number with 
eigenvalue a^ (E) t^ = —1. For generic smooth profiles 
of A, V and to, the index is a topological invariant that 
depends only on the asymptotic values of these functions. 
To see this we need only remember that the positive and 
negative energy states of h are paired: If 



then 



and 



htpE = Etpi 



h{a'^ (g) t'^Vb) = -E{a'^ t''Vb) 



1p-E =C7^ ® T^IpE 



If we deform the background fields in h smoothly, gener- 
ally the eigenvalues of h move and some modes that are 
at _B = could move away from E ~ and vice versa. 
However, since the non-zero modes come in positive and 
negative pairs, when modes move to or away from E = 0, 
they must do so in pairs. Furthermore, these pairs con- 
tain exactly one positive and one negative eigenstate of 



>T^: 



and 



{ipE + iP-e) = + (V'-B + iP-e) 



{ipE - iP-e) = - i'lpE - iP-e) 



Thus the difference between the number of zero modes 
with a^ ® T^ = +1 and a^ (E) t^ = —1, i.e. the index, is 
unchanged by a smooth deformation of the background 
fields. 

In the discussion above, we found that when either 
the mass term or the magnetic field Ba dominate the 
asymptotics, the index is 



Index(/i) 



(27) 



On the other hand, when the pseudo-magnetic field By 
is dominant at r — )■ c», the index diverges, as it should 
equal the total number of states in two infinitely degen- 
erate Landau levels. 



Index(/i) ~ - / (fxBv(x) 

TT J 



(28) 



(As we have discussed above, in order for Vg to be more 
important than the mass term at r — > oo, its fall off must 
be slow enough that the flux diverges at least logarithmi- 
cally in the volume.) This quantitiy diverges at least log- 
arithmically in the linear size of the 2-dimensional plane 
for the case in which it is valid and it would require more 
care to define it precisely. 

Another interesting topological quantity is the 7]- 
invariant. 



V{K)= Trsign(/i,) 



(29) 



where 



rh 



It is of interest because the state of a system of electrons 
which is governed by the single-particle Hamiltonian h^ 
and which has all negative energy levels occupied by elec- 
trons and all positive energ y le vels empty has electric 
charge given by the formula |17| 



partially cancel to leave a finite result. To make 77 un- 
ambiguous, we have added a term to the Hamiltonian, 
h ^ he ^ h + a^ (E) t^£- The term a^ (E) t^€ is physically 
equivalent to an out-of-plane component of the antiferro- 
magnetic order parameter which could be realistic. Since 



the square of /i^ obeys 



e^ >0 



and /le has no zero modes. All of its eigenvalues are 
either positive or negative. Furthermore, since no ma- 
trix anti-commutes with /i,;, the spectrum of h^ does not 
have particle-hole symmetry. Unpaired zero modes of h 
are also eigenfunctions of h^. They obey hipQ = and 
are eigenvectors of a^ E) t^ with eigenvalues -1-1 or —1. 
They are therefore eigenfunctions of h^ with energies e 
or — e. We can recover h by taking e to zero. Then, 
the spectrum becomes entirely symmetric except for un- 
paired zero modes, which were assigned to the positive 
or negative energy states oi h + a^ E> T^e. 

J] is a, topological invariant which is entirely determined 
by the asymptotic behavior oi A{x), V{x) and m{x). We 
can easily derive a formula for rj when we assume that 
the total flux 0y = -^ J (PxBv{x) is finite, but the to- 
tal electromagnetic flux cjiA = -^ J (PxBa{x) need not 
be. Computation of 77 is a straightforward generaliza- 
tion of that outlined in [3a] and also [la] , [13| and will be 
reviewed in the Appendix. The result is 



Vihe) = 



sign(e) 



Ve^Frr? 



2e 



V^^T: 



f^v (31) 



To derive this formula, we have had to assume that V{x) 
is such that the covariant derivative of the mass term, 
Vm(x) — iT^V{x)m{x) decays at least as fast as 1/r at 
large r. (For a vortex solution of Ginzburg-Landau ac- 
tion, this quantity generally decays exponentially.) 

Like the index, 77 in ([3T|) does not depend on the elec- 
tromagnetic gauge field A at all. There are a few inter- 
esting limits of ((3T|) . If we put e ^- 0+, we recover the 
index of h ([27]) . If we put By to zero, we see that 77 has 

which can bee attributed 



an irrational part, — , J . ,. 1 
to the polarization of continuum states by the vortex. 
Note that this polarization vanishes when e — > and the 
continuum spectrum becomes symmetric. The irrational 
part of the fractional charge of a vortex in the absence of 

MM- 



35 



V is identical to the one noted in Refs. [2£ 

Finally, if we decouple the vortex by setting the mass gap 

771 — >■ we recover 



(Q) = 7;Vihe) 



(30) 



77 is formally the difference of two infinite quantities, the 
number of states with positive energy and the number of 
states with negative energy. In interesting cases, these 



5. CONCLUSION 

We have discussed scenarios where vortices can have 
mid-gap electronic bound states. The main conclusion 



is that such bound states survive, practicahy unmodified 
in an external magnetic field. This in turn validates the 
consideration of vortices with mid-gap states as charge 
carriers in a Hall state. 

We comment that, like polyacetylene, where spin de- 
generacy obscures fractional charge, so that polyacety- 
lene actually has integer charged solitons, the system 
that we have discussed here has a valley degeneracy and 
the same phenomenon will occur. A vortex with a zero 
mode, actually has a pair of zero modes, one for each 
valley. Therefore, instead of being a two level system 
(zero mode occupied or unoccupied) where the states 
have charges -|-^ and — ^, it is a four- level system with 
charges —1,0,0,1. In polyacetylene, this leads to ex- 
otic spin-charge assignments, the charge ±1 states of the 
soliton are spinless, whereas the two charge states have 
spins ± 2 ' ^ topological version of spin-charge separation. 
A similar phenomenon would occur in graphene. When 
the mass condensate is a Kekule distortion, the vortex 
states would have spin degeneracy and the same exotic 
spin-charge assignments as occur in polyacetylene would 
now appear in two dimensions. If the condensate is an an- 
tifcrromagnct, the spin is replaced by the residual SU(2) 
valley quantum number, and could be more difficult to 
see experimentally, as the symmetry is emergent and the 
conservation law resulting from it is not exact. 

As for conductivity, as in polyacetylene, having vor- 
tices with low energy charged states makes the number 
of charged vortices available for conduction of electric 
current extremely sensitive to any displacement of the 
chemical potential from the charge neutral point. A sig- 
nature of charged vortices of this kind should be, like in 
polyacetylene, an enormous variation of carrier concen- 
tration with doping. 



which, using the fact that {cf^t^, h} — 0, we replace by 



2e ,. f°° , 
11 = — lim / dio Tr 



We define the current 



(7 T 



,-lih^ 



J{x) = tr {x\iGa^T^ ^°^ 



/l2 



'PK 



hl + e'^ + a;2 



\x) 



(A.2) 



(A.3) 



where tr indicates a trace over the Dirac and spin indices 
only. (Our previous trace, denoted Tr was a trace over 
Dirac and spin indices and position space as well TrO = 
j dx iv (x\0\x). ) Taking a divergence of the current 
yields the identity 



V • J = V • tr (xKggV ^, -\x) 



'0 



3,3 



tr {x\hQa r 



-2tr {x\a 



3^3, 



hae-f^'^i 



'«l^) 



3^3 



+2{e^ +w^)tr (a;|o-^T 



a^T^ 



.-UK 



hoe- ^''i ho 



\x) 



hl + e^+oj^' 



(A.4) 



Using the identity in ()A.4p . we rewrite (jA.2p as 
duj 2e 



rt = lim 

P^Ojo 



1 



Tr 



^3^3g-/3h^ 



3^3 



- / d^xV • tr {x\iaa'^T 



hoe-Pf'o 



:\x)\{A.5) 



The first term in (jA.5|) is evaluated using an asymptotic 
expansion of the heat kernel. 
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x) 



1 1 

47r/3 47r 

+& X Dm — m'^ 



[a^BA + r'^Bv) 

0{P) (A.I 



where Dm = (^ ^ 2iT^V\ m. When we multiply by 

a^T^ and take the trace, we see that the first term in 
(|A.5P is 



Appendix: Computation of the 77-invariant 



In this appendix, we outline the computation of ?/, the 
result of which was quoted in Eq. (|3T|) . We regulate the 
trace in ([29]) and represent the step function sign(/i) by 
an integral 



2 f°° 

ri = — lim / duj Tr 

IT p-^oJo 



,-pK 



h^ 



(A.l) 



n duj 2e 
lim / ir 



cjW-PK 



sign(e)- / d'^xB 



(A.7) 

Now, we consider the second term on the right-hand- 
side of (jA.Sp which, using Gauss' theorem, can be written 
as a surface integral on the circle at r = 00. To evaluate 
it, we must make some assumptions about the asymptotic 
behavior of the fields. Wc assume that the gauge field A 
is such that Ba grows no faster than a constant at large r. 
This allows the case of a constant magnetic field. In this 
case asymptotically, in a radially symmetric gauge Ag = 
-^r^. We shall need to assume that the gauge field V{r) 
decays at least as fast as 1/r for large r so that Dm ^ 1/r 



and By ^ t-- We then can use a functional Taylor 

expansion of the integrand in powers of Dm, where it 
turns out that only the first (linear) order survives the 
traces and the asymptotic and /? — >■ limits. The trace 
of the zeroth order vanishes. The first and higher order 
terms in the expansion are finite in the limit ^ — >■ which 
can now safely be taken. The first order in Dm is 



lim tT[x\ia<T T 



j2„, _3 



hi 



\x) 



tiidcr'W'^ / d y a'^m{x)g{x,y)a x Dm{y)g{y,x) (A. 8) 



where 



gix,y) = ix\- 



-D^ - (7^Ba - CT^T^B^ 



■UJ^ 



-\v) 



(A.9) 

Since the Green functions in (|A.8|) arc short ranged and 
|a;| — > oo, they will have support in the region where 
|y| — ?► oo and the background fields in (|A.9P can be re- 



placed by their asymptotic values. Then, the Green func- 
tion can be found explicitly, 



ff(^'y) "" ^ ' ^'^ 



-x(3:,y,A) 

sinhS^A 



with x(a;,y,A) = ^{^ - yfcoi\\BA\ + \(p'^Ba + 
jr? + w^ -I- e^) -I- i^x x y. We can now easily show 
that, because of the large r limit, the green func- 
tions ...g{x,y)...g{y,x)... in (|A.8P can be replaced by 
y). This can be seen by defining 



ix 



X = rx, rescaling y — )■ r y and taking the large r limit. 
Then, we have 



lim trice kr • aa t —^ -\x) 

1 e^^m^dem^ + Im^Ve 



ni^ + o;^ -f- e^ 



(A.IO) 



Doing the remaining integrals, wc obtain pip . 
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